Generation of pure single photon wavepackets by conditional preparation based on 

spontaneous parametric downconversion 



O 
O 

(N 

> 
O 

(N 



Alfred B. U'Ren 1,2 , Christine Silberhorn 1 , Reinhard Erdmann 3 ,Konrad Banaszek 1 , 
Warren P. Grice 4 ,Ian A. Walmsley 1 and Michael G. Raymer 5 
Clarendon Laboratory, Oxford University, Oxford, 0X1 3PU, UK 
2 Division de Fisica Aplicada, Centro de Investigacion Cienttfica y 
Educacion Superior de Ensenada ( CICESE), Baja California, 22860, Mexico 
3 Rome Air Force Laboratory, Rome, NY, USA 
4 Center for Engineering Science Advanced Research, Computer Science and Mathematics Division, 

Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
5 Department of Physics and Oregon Center for Optics, University of Oregon, Oregon 97403, USA 

(Dated: February 1, 2008) 

We study the conditional-preparation of single photons based on parametric downconversion, 
where the detection of one photon from a given pair heralds the existence of a single photon in 
the conjugate mode. We derive conditions on the modal characteristics of the photon pairs which 
ensure that the conditionally prepared single photons are quantum mechanically pure. We propose 
specific experimental techniques which yield photon pairs ideally suited for single photon conditional 
preparation. 
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I. INTRODUCTION 

Pure single photon states are probably the most fun- 
damental entities in quantum optics, and constitute the 
starting point for many optically-based quantum en- 
hanced technologies. A basic requirement for many key 
applications is the ability to generate reliably pure sin- 
gle photon wavepackets capable of high-visibility inter- 
ference. Single photon wavepackets may be generated 
from number-correlated pairs by conditional state prepa- 
ration. For example, photon pairs produced by the pro- 
cess of spontaneous parametric downconversion (PDC) in 
a v/ 2 ) non-linear medium allow conditional preparation 
of single photons, since the detection of one photon in the 
pair heralds the presence of the conjugate photon. PDC- 
generated photons have been employed for both funda- 
mental tests of quantum mechanics and to demonstrate 
various quantum communication applications, such as 
quantum teleportation[l|, quantum dense coding [2| and 
quantum cryptography [3j- 

A recent proposal for linear optical quantum compu- 
tation (LOQC)[4| exploits the photon bunching that oc- 
currs in quantum interference between pure state single 
photons in conjunction with conditional state prepara- 
tion in linear optical networks. There the detection of 
auxiliary photons indicates the successful operation of 
a given gate. For LOQC, as for all other applications 
involving optical quantum networks, high-visibility in- 
terference between photons from multiple sources is cru- 
cial. This necessitates, on the one hand, precise timing 
of the photons. This, in turn, demands a high proba- 
bility of simultaneous generation from synchronized dis- 
tinct sources. On the other hand, precise control of the 
modal structure of the generated photon pairs is essen- 
tial to guarantee the purity of the conditionally prepared 
photon states. For PDC based implementations these re- 



quirements lead to the need for a pulsed pump and for 
efficient sources exhibiting high brightness and detection 
efficiencies. Spatial modal control of the emitted pho- 
tons, and a much increased production rate of photon 
pairs may result from waveguided PDC[5|. 

A common approach to guarantee indistinguishability 
in optical experiments uses strong spatial and spectral 
filtering. The cost of this, however, is a prohibitive re- 
duction of the generation rate of usable single photon 
wavepackets. Since common birefringent phase-matching 
leads to spatial-spectral correlations with overlapping 
contributions, spatial filtering also causes optical losses 
for the signal photon for finite spectral filter bandwidths. 
As a consequence the fidelity of the conditionally pre- 
pared signal photons is limited due vacuum contribu- 
tions. To eliminate these and to ensure, for example, 
successful gate operations in LOQC post-selecting by co- 
incidence measurements is inevitable. But such post- 
selection, in turn, hinders the realization of scalable net- 
works involving more than one or possibly two sources. 
Waveguided PDC is one route to accomplish both, to 
decorrelate the spatial degree of freedoms from spectral 
ones and to eliminate at the same time the spatial corre- 
lations between the PDC photons. This implies no need 
for spatial filtering, however, conditionally prepared sin- 
gle photons are in general still described by a mixed state 
due to spectral correlations present in the photon pairs. 
In this paper we study the conditional state preparation 
of single photon wavepackets based on PDC pumped by 
ultrafast pulses. We use a multi-mode description and 
propose novel engineering techniques for eliminating cor- 
relations between unobserved variables. This ensures 
that the detection of one photon from a given pair di- 
rectly yields a pure single photon state in the conjugate 
mode, without resorting to filtering. The latter consti- 
tutes an important step for the further development of 
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quantum networks since it is expected to lead to an in- 
crease in the probability of simultaneous generation from 
multiple crystals by several orders of magnitude. 



II. CONDITIONAL STATE PREPARATION OF 
PURE SINGLE PHOTON STATES IN A 
MULTIMODE DESCRIPTION 

Single photon states are defined as containing a sin- 
gle quantum in the photon number basis. A complete 
specification of the photons requires that all degrees of 
freedom associated with the quantized electromagnetic 
field of the occupied optical modes must be taken into 
account. Thus, the complete description of the photons 
involves their polarization /i, their spatial vector r or 
the wavevector k, and their frequency u> or time depen- 
dence t. The dispersion relationship between k and u> 
reduces one degree of freedom, such that we can restrict 
our analysis to transverse wavevectors k^. For simplic- 
ity we also consider only one single polarization, which is 
appropriate for PDC experiments aimed at conditional 
prepararation of single photons. A pure single photon 
wavepacket then corresponds to a coherent superposition 
of monochromatic plane waves: 



= N 



J dui J dk c(uj,k) exp[i(k-r — cut)] ^(ujjk) |0), 

(1) 

where N denotes the normalization constant and the co- 
efficient function c(u>, k) determines the modal structure 
of the photon. The corresponding density operator is 
given by p = \(f>)(4>\ and exhibits nonzero off diagonal 
elements. This is to be contrasted with a completely 
mixed single photon state in frequency and momentum 
for which no coherence between the different modes exists 
and is characterized by a density operator of the form 

p = N J doj J dk\c{u,k)\ 2 a\oj,k) |0> <0j a(u),k), (2) 

which does not exhibit any off-diagonal elements. 

In the process of PDC individual pump photons split 
into two daughter photons in modes a and b. For the 
purpose of conditional state preparation the modes d 
and b should be distinguishable in at least one degree 
of freedom, which might be polarization for type II 
phase-matching or the spatial degree of freedom for non- 
collinear phasematching. The interaction Hamiltonian of 
the x^-process for a strong pump field, which is treated 
classically, results in state of the approximate form: 

|*) = |0> |0> + n J J dw.dw, J J dk^dkt 

xFK,^;^,^) at(c s ,k s ^) b^k^) |0) (3) 

where fi is a measure of the efficiency of the PDC 
process and the function F(u s , k^-; w i; k^~) represents a 
weighting function for the different spatial and spec- 
tral modes present, which results from the pump enve- 
lope and phasematching functions defined by the spe- 
cific downconversion configuration [8]. It is assumed that 



states of higher photon numbers can be ignaored. The 
function F(cu s , k^~; u>i, k^) may contain correlations be- 
tween the signal and idler photons in the spatio-temporal 
continuous degrees of freedom, thus yielding an entangled 
two-photon pair. A useful tool for analyzing such entan- 
glement is that of carrying out a Schmidt decomposition 
[9j of the spatio-temporal state amplitude. The Schmidt 
decomposition entails expressing the two photon state in 
terms of complete basis sets of orthonormal states AjjO) 
and B^\0) such that: 



l*) = Ev / ^4 t i ^|o) 



(4) 



where X n — 1 and where the values A„ are known 
as the Schmidt coefficients. In this description the ef- 
fective signal and idler creation operators A^, B^ are 
given in terms of the Schmidt functions ip n (u> s ,kjr) and 
(f> n (ui,kl) as 



dki 



and 



Bl = / dk 



dw s M^ s ^s) a f (w s ,kj-)|0) (5) 



while the two-photon state amplitude F(kjr, u> s ; k^, oji) 
can be expressed as 

F(uj s ,k^;uj u kl) =^^4(^,^)^(01,,^), (7) 



where u a i and k^ represent the wavelength and the 
transverse wavevector of the signal and idler photons. 
For the calculation of the Schmidt functions themselves 
the first step is to determine the reduced density matrices 
of the two sub-systems corresponding to the signal and 
idler photons, which are defined by 

Ps (LU s ,ki;Cj s ,ki)= fjdoj' dk ± 'F(oj s ,kj,u; , ,k ± ') 

xF*(^,ki,cy,k^') (8) 

and 

p^ki-Co^ki) =JJduj' dk J - / F(w / ,k J - / ,w i ,k i L ) 

xF^'.k 1 ,^,^') (9) 

and as a second step solve the following the following 
integral eigen- value equations Q: 

JJ duj'dk^'p s (co s , k^, «/, k^')^ n (co', k^') 

= A„V>„(w s ,k^) (10) 

// dw'dk^piiuji, kf,w', k x/ )<£ n (u/, k x ') 

= \Mio ll ki) (11) 
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thus yielding the Schmidt functions ipn and 4> n and the 
joint eigenvalues A„ as Schmidt coefficients. The Schmidt 
decomposition simplifies the quantum state description, 
by turning the original state expressed as a quadruple 
integral over frequency and transverse momentum vari- 
ables into a discrete sum, typically containing only a lim- 
ited number of terms. The Schmidt functions ip n (u> s , kjr) 
and 4> n (uii,]i^-) can be thought of as the basic building 
blocks of entanglement in the sense that if the signal 
photon is determined to be described by a function ip n , 
we know with certainty that its idler sibling is described 
by the corresponding function </>„. The Schmidt modes 
also represent space-time localized wavepackets [6|, and 
may be thought of as wavefunctions of a photon Q. 
The probability of extracting a specific pair of modes 
[^ n (w s ,k^),0„(wi,k^)] is given by the parameter A„, 
which is real and nonnegative. The amount of entangle- 
ment can be conveniently quantified by the cooperativity 
parameter (To| defined in terms of the Schmidt eigenval- 
ues as: 



K = 



J2 



(12) 



The value of K gives an indication of the number of ac- 
tive Schmidt mode pairs, which in turn is a measure of 
how much entanglement is present in the photon pairs. A 
two-photon state for which the cooperativity assumes its 
minimum allowed value K — 1 represents a state in which 
there is a single pair of active Schmidt mode functions 
and therefore exhibits no spectral or spatial entangle- 
ment. Another measure of the degree of non-separability 
is the so-called entropy of entanglement: 



s = -$> n log(A n ) 



(13) 



which vanishes for a factorizable state, and increases 
monotonically with the amount of entanglement present. 

For the conditional state preparation of single photon 
states via parametric down-conversion the idler photon is 
treated as a trigger, such that in the ideal case of perfect 
photon number correlations, a detection event of the trig- 
ger idler channel heralds the emission of a photon in the 
signal channel. The quantum state resulting from this 
conditional preparation can be obtained by modelling the 
trigger detection process in terms of a projection opera- 
tor and tracing over the trigger channel photon: 



p s = Tr^pfi). 



(14) 



Here the subscript i denotes a partial trace acting over 
the trigger idler mode, p is the density operator describ- 
ing the two-photon PDC state, while II is the measure- 
ment operator modelling the trigger photon detection. 
From the Schmidt decomposition it can be seen eas- 
ily that filtering one specific Schmidt mode </>„ (uii , k x ) 
with the projection operator B^\0)(0\B n allows condi- 
tional preparation of a pure photonic wavepacket in the 



signal mode. However this procedure requires coherent 
frequency-time and spatial filters, which are non-trivial 
to implement. 

We thus proceed to investigate the synthesis of pure 
signal photons by time integrated detection of the trig- 
ger idler photon in combination with passive spectral and 
spatial pinhole filters. Correlations between the signal 
and idler photons generally lead to incoherence between 
different frequency and spatial components upon detec- 
tion, which gives a mixed output signal state. The phys- 
ical reason for this is that a measurement of the trigger 
reveals some information about the properties of the idler 
and hence destroys indistinguishability. A spatial and 
time-integrated detection of the idler state with a spec- 
tral interference filter with a profile function <j(u>) and a 
spatial filter of £(fc) can be modelled by the projection 
operator 

II = J dt J dr J J duj dd> J J dk dk^a{uj) £(k x ) 

x exp[i(k • r - wt)]ft(w, k) |0)i(0|;S(cD, k) 
x expHQc • f - ut)] cr*(cD)£*(k x ) 

= / 'dcu J dk^ \a{u>)\ 2 \^)\ 2 &t( w , k ) |0>«<0U 
x6(w,k). (15) 

By carrying out the calculation expressed by Eq. ITU using 
the density operator |^ r )(4 r | of the PDC state, it can be 
shown that the quantum state describing the signal chan- 
nel, upon registering a trigger detection event, is given 

by 



da; / dk- 1 / dw, / dk^ / du>, / dkj- 



Ps -- 

x|aH| 2 |e(k x )| 2 F(uj s ,k^;uj,k x ) tf(u.,\£) 
x\0) s (0\ s a(u> s ,k^)F*(uj s ,k^-,Lj,k x ) . (16) 

If we utilize the Schmidt decomposition of Eq. \7\ the den- 
sity operator of the conditionally prepared signal photon 
can be cquivalcntly rewritten as 

p s = J^V^A^ [du> [ dk 1 - \a(Lj)\ 2 le^ 2 

n m J J 

x^k^C^k^itjO^oU™ (17) 

Under what conditions does p s represent a pure state? 
Purity requires that the density operator p s contain a 
single term, such that the signal state can be expressed 
in terms of a single Schmidt mode pair. From the 
expressions in Eq. [16] and Eq. [TT] two approaches to 
achieve purity for the conditional state preparation via 
PDC become apparent. Restricting the trigger detec- 
tion bandwidth to ideally one single spectral component, 
i.e. |<t(w)| — * 5{oj — fi), as well as applying strong spa- 
tial filtering with an imaging lenses and a narrow pin- 
hole, i. e. |£(k )| — > 8(k — K) effectively suppresses 
the integration over u> and k^ in Eq. [TB] and eliminates 
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the incoherent sum over different frequency and spatial 
components, leaving the signal photon in a monochro- 
matic plane wave. Note, however, that this method for 
purification of the prepared signal output state necessar- 
ily implies diminishing count rates and approaches truly 
pure states only in the limit of vanishing counts. The 
Schmidt decompositions suggests an alternative route for 
the preparation of pure signal photon states exploiting 
PDC pair generation. If there exists only one Schmidt 
pair present, the sums in Eq.[T7]disappear and the double 
integral term becomes an overall multiplicative constant. 
The latter implies that the prepared photon now forms 
a truly quantum-mechanically pure state. For a better 
understanding of this result we can relate the Schmidt 
decomposition to the correlations between the generated 
photons. The existence of only one Schmidt mode is 
actually equivalent to the statement that the PDC two- 
photon state is factorizable. Under such circumstances, 
the spectral and spatial properties of signal and idler pho- 
ton are independent from each other so that the detec- 
tion of the idler photon yields no information whatsoever 
about the signal photon. In order to further quantify the 
purity of the conditionally prepared single photons in the 
signal arm we can introduce the purity parameter p de- 
fined as 

p = Tr(p2). (18) 

Evaluating the purity p in terms of the Schmidt decom- 
position, we can clarify the connection between purity 
of conditionally prepared single photons and the entan- 
glement present in PDC photon pairs. In the absence of 
filters p = J2 n ' which is equivalent to p — -i , where 
the cooperativity parameter K quantifies the entangle- 
ment of PDC. Thus, entanglement between the signal 
and idler photon hinders pure state preparation or, con- 
versely, spatial and spectral separability (which can be 
enhanced by appropriate filters) ensures purity for con- 
ditional state preparation setups. This criteria delivers a 
valuable tool for engineering of PDC photon pair sources 
optimized for the generation of pure single photon states 
with without resorting to strong filtering and illustrates 
the necessity of decoupling the signal and idler photons 
in all degrees of freedom. 

III. SPECTRALLY DECORRELATED 
TWO-PHOTON STATES 

As shown in section UIl the generation of pure single 
photon wavepackets based on PDC photon pairs neces- 
sitates that all correlations between the signal and idler 
photons be eliminated. This includes correlations in all 
degrees of freedom including spectral, transverse momen- 
tum and polarization. Since the full discussion of all de- 
grees of freedom is beyond the scope of this paper, for the 
discussion following we assume that the spatial modes of 
the PDC photon pairs are independently rendered decor- 
related. Such spatial decorrelation can be achieved, for 



example, by waveguided PDC[f|, in which the photons 
may be emitted only into specific transversely-confined 
modes. In addition, some of the techniques to be pre- 
sented here designed to obtain spectral decorrelation for 
bulk crystals, could be extended to the spatial domain. 

A number of techniques have been proposed to gener- 
ate photon pairs exhibiting spectral decorrelation with- 
out resorting to spectral filtering. Grice et al. showed flT| 
that an avenue towards such states is the group velocity 
matching condition derived by Keller et a/. [121]. Such 
a technique is extended here, showing that if an addi- 
tional condition between the pump chirp and the crys- 
tal dispersion is fulfilled, factorizability of the joint spec- 
tral amplitude is guaranteed. Considering the complete 
joint amplitudes rather than restricting attention to in- 
tensity correlations is important, because correlations in 
the phase terms of the joint spectral amplitude, which are 
not apparent in the joint spectral intensity, can introduce 
correlations in the times of emission. Thus, in order to 
guarantee full signal-idler decorrelation, the joint spec- 
tral amplitude (and therefore the joint spectral intensity) 
must be factorizable. 

It has likewise been shown that a spectrally decorre- 
lated state may be obtained by exploiting transverse mo- 
mentum in a bulk x' 2 ' type-I crystal[lj|. Specifically, 
this is achieved using of a focused Gaussian pump beam 
with a spot size fulfilling a certain relationship with the 
crystal length. Walton et al. have recently reported 
a different scheme in which a transverse pump gener- 
ates counter-propagating spectrally de-correlated photon 
pairs [II]. The essential advantage of these two tech- 
niques is that spectral decorrelation may be obtained at 
any wavelength where phasematching is possible, as op- 
posed to the technique based on group velocity matching 
which occurs for quite restrictive frequency regimes. 

In this paper we also introduce a novel technique, 
based on a sequence of crystals with intermediate bi- 
rcfringent spacers. In this scheme — although group ve- 
locities are not matched to each other — the maximum 
group velocity mismatch observed can be limited to an 
arbitrarily small value. Furthermore, spectral decorrela- 
tion, as well as a more general class of spectrally engi- 
neered photon pairs, may be obtained at any wavelength 
where phasematching is possible [lj|. 



A. Group velocity matching 

The general expression for the joint probability ampli- 
tude of the two-photon state generated in the process of 
PDC in the ultrashort pulsed pump regime (and assum- 
ing fixed directions of propagation^ given by 

f(uj s , uji) = a(uj s ,uJi)cj>{uj s ,u)i) (19) 

where the pump envelope function a(w s ,o>.;), modelled 
here by a Gaussian in terms of the frequency detunings 
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from the central PDC frequency Un = w„ — ujq: 



a(v s + Vi) = cxp 



V s + Vi 



(20) 



and where the so-called phasematching function (j)(u> 8 ,uji) 
describing the optical properties of the nonlinear crystal 
is expressed as[T 



' LAk((jj s ,oJi)\ jiAi 
{u) s ,oji) = sine I I e 



(21) 

in terms of the phase mismatch: 

Ak(u; s ,LUi) = k s (uj s ) - ki(un) - k p {uj p ). (22) 

We will at this point restrict our attention to a 
frequency-degenerate process, i.e. the Taylor expansion 
for each of the two photons is assumed to be centered at 
the same frequency, luq. We can then express the-phase- 
mismatch Eq. [22] as a Taylor expansion up to second or- 
der: 

LAk(v s , Vi) = LAk a + t s v s + nvi 

+p s v 2 s + fcvf + (3 p v s v l + 0(v 3 ) (23) 



with 



Afc° = k s (uj Q ) + ki(u ) - k p (2uj ) 



(24) 



representing the constant term of the Taylor expansion, 
which must vanish to guarantee phase-matching; 0(v 3 ) 
denotes terms of third and higher order and (// = s,i) 



t^ = L [fc> ) - *£(2wo)] = 
/3 P = Lk;(2Lu ). 



(25) 



In the previous expressions ' and " denote first and second 
derivatives with respect to frequency(evaluated at lu in 
the case of the signal and idler wave vectors and at 2wo 
in the case of the pump wavevector). u p represents the 
pump group velocity whereas it M (/i = s,i) represents the 
group velocity experienced by each of the signal and idler 
photons. 

We will use the additional approximation of express- 
ing the sine function in the phasematching function as a 
Gaussian function through the approximation: 



sine (a;) 



with 7 = .193. 



(26) 



where the numerical value of 7 results from the condition 
that both functions exhibit the same FWHM. We will 
further allow the pump field to carry a quadratic phase 
to account for chirp. Under such circumstances, the joint 



spectral amplitude is given as: 

f(v s > Vi) = 

M exp 
x exp 
x exp 



exp [ij3 t {v s + ViY 



--^{TsVs + TiVi) 2 



(27) 



i-^(T s v s + nVi + [3 s v 2 s + (3 t vf + (3 p v s Vi) 



where M is a normalization constant and (3t is the GVD 
dispersion term which the pump experiences prior to the 
crystal. We note that the above expression contains, 
within the regime of the approximations used all terms 
up to quadratic order in the frequency detunings. By ex- 
panding the exponential terms we can express the joint 
spectral amplitude as 



f(v Sl Vi) oc 



exp 
x exp 
x exp 



7 2 \ 2 



7 2 \ 2 



exp 
exp 



1-^V S + l[ p t + —V 



■ Ti , ■ la 1 & 2 
i-^Vi + i[Pt + —v, 



2 

2 



-2 I -4 + T - -") ''-''< + ' ( 2/3* + -r- I i',v, 



.(28) 



Here the contributions depending on v s and Vi exclusively 
and those depending on both v s and Vi are written as 
separate factors, and phase terms linear in v^ have been 
dropped. We can see from Eq. [55] that any unfactoriz- 
ability resides in the mixed term containing an argument 
proportional to v s Vi. Furthermore, we can now easily see 
that the conditions that guarantee a factorizable state 
(i.e. that make the mixed term vanish) are: 



7T s Tj = 



and 



2(3 t + 



0. 



(29) 



(30) 



We note that the second condition [Eq. [30] refers to 
phase contributions in the joint spectral amplitude, and 
was not included in the earlier work by Grice et aZ.[Tl|. 
Because such phase contributions play no role in the in- 
tensity, the fulfilment of the first of the derived conditions 
[Eq. [55] is sufficient to obtain a factorizable joint spec- 
tral intensity. If both conditions [Eqns. [5pJ and [3D] are 
satisfied, the two-photon state is guaranteed to be free 
of correlations in the spectral and temporal domains, as 
indicated by a factorizable joint spectral amplitude. We 
now analyze the effects of the first condition on the two 
photon state. If condition Eq. [55J is fulfilled, we may 
express the joint spectral intensity as: 



(31) 



S{uj Sl uji) = \ f(uj Sl uJi)\ cx exp 









oc exp 




exp 
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where the spectral widths a s and Ci are given by: 

2(7 



4 + 7<t 2 t2 



(32) 



with fj, = s,i. Thus we have confirmed that such a joint 
spectral intensity exhibits the desired factorizability. The 
ratio of the larger of the two spectral widths (for signal 
and idler) to the smaller one, to be referred to as the 
aspect ratio, (i.e. a measure of the degree of elonga- 
tion exhibited by the two-photon spectral distribution in 
{u s ,u>i}) is given by: 



7(7 2 T 2 2 ' 



(33) 



where t\ = t s and r 2 = if a s > Ui and likewise t\ = t% 
and T2 = t s if ai > a s . In what follows the aspect ratio 
will prove to be a useful tool in the characterization of 
the two-photon state. 

Let us now consider the effect of the second of the 
derived conditions [Eq. [3D] . We will start by analyzing 
the effect of the phase terms present in the joint spec- 
tral amplitude [Eq. 128] on the photon pair correlations. 
First, the linear phases ir^v^ (with /j, = s,i) shift the 
times of emission with respect to the pump pulse for 
each of the two single photon wavepackets. The fact 
that for a type-II interaction t s ^ Ti implies that the 
two single-photon wavepackets forming a given photon 
pair are emitted at different mean times, where the mean 
emission time difference is simply r s — Tj . This temporal 
walkoff can be compensated by means of a relative delay 
between the signal and idler photons introduced after the 
crystal. While the GVD terms (/3 S ,A and fi p ) have no 
effect on the joint spectral intensity of the two-photon 
state (in the regime of the approximations used which 
neglects contributions of cubic and higher order terms in 
the frequency detunings) they translate into broadening 
in the temporal domain, which becomes apparent upon 
a Fourier transformation of the joint spectral amplitude 
to obtain the joint temporal amplitude. The quadratic 
phase i(Pt + Pn/2) v %, (with fx = s,i) results in uncorre- 
cted temporal broadening of the two photon wavepackets 
along the t s and ti axes. The temporal broadening is pro- 
portional to the crystal length and is, in general, differ- 
ent for each of the two photons [20j (since (3 S ^ /%). The 
mixed term, proportional to u s Ui, introduces correlations 
in the times of emission. This means that decorrelations 
of spectral intensity are not sufficient for true single mode 
photon pair generation. We note that it is the (3 P term 
which is responsible for extirpating the factorizability. 
Therefore, if it were possible to eliminate the effect of /3 P , 
the two-photon joint temporal amplitude would broaden 
in such a way that the decorrelated character is main- 
tained. Fortunately, it is possible to compensate for the 
presence of (3 P simply by letting the incoming pump field 
have a chirp (quadratic phase) which fulfils our second 
derived factorizability condition (Eq. |3"U|) . The quadratic 



phase which the pump should carry is therefore given by: 
Pt = — /3 p /4. We also note that it is in principle possible 
to "tune" the degree of temporal correlation in the two- 
photon state (and therefore the level of distinguishability 
in an interference experiment) by using a pump field with 
a variable chirp, i. e. variable /3f 

Let us return to the first of the conditions derived for 
factorizability [Eq. I29| . This condition may be written in 
terms of the crystal and pump field parameters as: 



+ 1 L 2 {k' s -k')(k' l -k')=0 



(34) 



where a is the pump bandwidth, L is the crystal length, 
k' p represents the first frequency derivative of k p evalu- 
ated at 2wo, (/i = s,i) represents the first frequency 
derivative of fc M evaluated at loq and 7 « .193. For spe- 
cific experimental situations, pairs of values for the pump 
bandwidth a and the crystal length L may exist such 
that the above condition (Eq. IM]) is fulfilled. Note that 
for the condition in Eq. [34] to be fulfilled, one of the fol- 
lowing must be true: k' s < k' < k[ or k[ < k' p < k' s , i. e. 
the group velocity of the pump must lie between that of 
the signal and idler. None of the experiments recently re- 
ported in the literature make use of crystals meeting the 
above requirements. This should not be surprising, since 
it requires a material in which one of the daughter pho- 
tons has a smaller group velocity than the pump, which 
necessarily has a shorter wavelength. Alternatively, this 
condition may be expressed as the requirement that t s 
and Ti have opposite signs. 

In order to get further physical insight, consider the 
phasematching condition resulting from the approxima- 
tion of neglecting second and higher order terms in the 
Taylor expansion (and assuming that the constant term 
vanishes). In this case the phasematching condition [see 
Eq. [22 becomes: 



TiVi = 



(35) 



The contour in {v s -Vi\ space defined by perfect phase- 
matching is therefore a straight line with slope —T s /ri. 
In other words the angle subtended by the v s axis and 
the perfect phasematching contour is given by: 



'II 



-arctan 



= —arctan 



k e 



(36) 



where the II subscript refers to type-II phasematching. 
Note that the angle 9jj does not depend on the crystal 
length; it is solely a property of the dispersion of the 
material at the particular PDC wavlength. For type-I 
PDC, t s = Ti so that a similarly defined angle has the 
fixed value 9] = —45°. Thus, the condition that k' s — k' p 
and k\ — k' p have opposite signs cannot be met with type- 
I PDC while in the case of type-II it translates into the 
requirement that the slope of the perfect phasematching 
contour be positive. 

Let us now consider a specific case fulfilling the positive 
slope requirement: the case where the slope of the perfect 
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phasematching contour is unity, i.e. Qpu = 45°. It is 
straightforward to show from Eq. [35] that the condition 
which guarantees such a unit slope is given by: 

V -+- V 

^p = fc; (37) 

which says that the average inverse group velocities for 
the signal and idler photons equals the pump inverse 
group velocity. This condition, which can be alterna- 
tively expressed as t s + Ti = OAs referred to as the group 
velocity matching condition [ill, E3] • We can picture this 
condition as the requirement that the signal and idler 
photons are temporally delayed with respect to the pump 
in a symmetric way, with the two PDC photons delayed 
in opposite directions, by the same amount, from the 
pump. Alternatively we may interpret the group velocity 
matching condition Eq. [37] for a factorizable two-photon 
state in terms of the joint spectral intensities and identify 
that the aspect ratio [see Eq. [33] evaluates to unity. The 
latter implies that the factorizable joint spectral inten- 
sity exhibits a circular shape on the v a — plane. Such 
symmetry is crucial for experiments which rely on clas- 
sical interference between signal and idler PDC modes, 
since it can then fulfil the requirement of ideal signal-idler 
spectral mode matching and lead to unit visibility in the 
conventional Hong-Ou-Mandel experiment [Hj]. Hence in 
this configuration the generated PDC two photon state 
exhibits signal-idler mode matching as well as the factor- 
izability required for ideal conditional preparation of pure 
single photons, and therefore constitutes a versatile two- 
photon state, useful for a broad class of experiments [3] • 

B. Asymmetric group velocity matching 

In the context exclusively of generating single pure sin- 
gle photon states, each of the crystals used must generate 
factorizable (but not necessarily symmetric) two-photon 
states in order to guarantee unit interference visibility 
in a multi-crystal experiment where only signal or re- 
spectively idler photons are superimposed on each other. 
Thus, in this context, a non-unit aspect ratio is indeed 
not a limitation (as long as the elongation is "aligned" 
with the frequency and time axes). In this section we 
study experimental situations where it is in fact desir- 
able to obtain a state exhibiting a spectrally "elongated" 
shape and thus a high aspect ratio. 

The first condition for frequency decorrelation (Eq. [29]) 
can in fact be met for several common \^ crystals. 
For example, in the case of periodically-poled KTP (PP- 
KTP) (in degenerate collinear operation), it can be met 
for configurations where the PDC central wavelength is 
within the range 1.207/im < A < 2.364/im (see also 
Ref. [IB]). Within this range only the particular wave- 
length of A = 1.568/mi fulfills the group velocity match- 
ing condition (Eq. I37|). which yields the unit-aspect ratio 
two-photon state, while the aspect ratio departs from 
unity for all other wavelengths. The lower and upper 



bounds correspond to two-photon states characterized by 
a high aspect ratio with asymmetric joint spectral inten- 
sities, though maintaining "alignment" with the v s and 
Vi axes). In Ref. [ll[, it was similarly established that in 
the case of BBO factorizability may be obtained for wave- 
lengths between 1.169nm and 1.949nm with the group ve- 
locity matching condition fulfilled at 1.514nm. This spec- 
tral regio is of great importance for optical quantum com- 
munication using optical fibers. Unfortunately for the 
materials discussed above, frequency-decorrelated PDC 
(and specifically unit aspect ratio decorrelated states) oc- 
curs at wavelengths at which single-photon detectors are 
not well-developed. Although InGaAs avalanche photo- 
diodes have been shown to be capable of photon counting 
at such longer wavelengths fl7j . quantum efficiencies are 
rather low (about 15%) while dark count control neces- 
sitates cryostatic cooling. Thus, frequency decorrelated 
two-photon states, which occur for available crystal ma- 
terials at longer > 1/xm wavelengths, present acute exper- 
imental challenges mainly in terms of detection. In what 
follows we will discuss the synthesis of factorizable two- 
photon states in arbitrary spectral region, with a view to 
accessing an experimentally more convenient range. 

We begin our analysis by re-writing the condition for 
factorizability (see Eq. I3~i|) as: 

— + 7(77-, = 0. (38) 

Consider this condition in the long PDC crystal regime. 
Recall that the temporal walkoff terms r s and n are each 
proportional to the crystal length L. Thus, we see that if 
t s >> (as is the case in the limit L — ► oo), the con- 
dition reduces to the simpler constraint: t, : = 0: i.e. the 
latter tells us that spectral decorrelation can be achieved 
employing a long crystal while making one of the tempo- 
ral walkoff terms vanish. Because in this variant of the 
technique the pump group velocity is matched to that 
of one of the PDC photons (but not to both) we refer 
to this technique as asymmetric group-velocity match- 
ing. Note that for a given crystal length L the frequency 
decorrelation condition [Eq. 134] can in general be met 
at most for one particular value of the pump bandwidth 
<r; i.e. there is a strict one to one relationship between 
crystal length L and the pump bandwidth a required for 
spectral factorizability. The latter can be an experimen- 
tal limitation if for example the pump bandwidth cannot 
be easily modified. In the long crystal regime, however, 
the condition becomes independent of the pump band- 
width a. The fact that a frequency decorrelated two- 
photon state may thus be synthesized irrespective of the 
pump bandwidth makes such a long crystal PDC a more 
flexible source of spectrally-engineered photon pairs. In 
such a long-crystal source, varying the pump bandwidth 
merely alters the aspect ratio without affecting the spec- 
tral decorrelation. 

Fortunately it is possible to obtain a decorrelated two- 
photon state (albeit with a high aspect ratio) at wave- 
lengths suitable for room-temperature operated silicon- 
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FIG. 1: High aspect ratio decorrelated two-photon state ob- 
tained by degenerate collinear PDC from a 2cm-long KDP 
crystal exhibiting a cut angle of 67.77°. (A) "Vertical" 
phasematching function exhibiting asymmetric group veloc- 
ity matching (B) Pump envelope function describing a pump 
pulse centered at 415nm with a spectral bandwidth of 5nm. 
(C)Resulting factorizable two-photon state. 



based on shorter crystals. 5) KDP crystals at a ~ 70° 
cut angle exhibit a relatively small transeverse walkoff; 
at a pump wavelength of 415nm the walkoff angle is 1.15° 
giving, for a 2cm long crystal a lateral displacement of 
just over 400//m. We note, however, that used as a con- 
ditional source of single photons in networks where signal 
and idler modes are not superimposed, the walkoff exhib- 
ited is not a limitation, as the signal and idler photons 
need not be indistinguishable from each other. 6) As will 
be discussed in the next section, such a high aspect ra- 
tio factorizable state additionally exhibits an interesting 
dispersion insensitivity effect. 

C. Dispersion insensitivity in long-crystal 
asymmetrically group velocity matched photon pairs 



based avalanche photodiodes. These detectors can ex- 
hibit both ~ 60% quantum efficiency as having a man- 
ageable level of dark counts. Figure[T]shows a two-photon 
state produced by collinear degenerate type-II PDC in a 
2cm-long KDP crystal cut at a phasematching angle of 
68° yielding photon pairs centered at 830nm. Note that 
although the spectral distribution is elongated, factoriz- 
ability is exhibited and hence pure conditional pure state 
generation (see Sec. [TTJ) is ensured. It is also interesting 
to notice that — in contrast to general PDC — such a state 
is predicted to yield conditionally prepared photon states 
which are Fourier transform limited The two-photon 
state features a narrow spectrum for the ordinary-ray 
photon and a broader one for the extraordinary-ray pho- 
ton. This arises from the matching of the group velocity 
on the pump photons to that of the o-ray photon, but 
not to the e-ray photon. 

In summary, our proposed long-crystal KDP-source 
yielding spectrally factorizable photon pairs, based on 
asymmetric group velocity matching, is characterized by 
the following properties: 1) If it is used for conditional 
single photon generation the prepared photons are de- 
scribed by quantum mechanically pure states. This prop- 
erty is evident in the by Fourier transform limited char- 
acter of the photonic wavepackets. 2) Narrow spectral 
bandwidth ordinary-ray photons allow spectral filtering 
to eliminate any background light without discarding 
PDC photons. If the filtered mode is then used as a 
trigger in a conditional source of single photons, the re- 
sulting photon exhibits the largest spectral bandwidth 
allowed by the pump bandwidth. 3) The decorrelated 
character is independent of the pump bandwidth, thus 
eliminating an important experimental constraint [21| . 4) 
The generated state requires a pump centered at 415nm, 
easily accessible with a frequency-doubled Titanium sap- 
phire laser, while the produced photons at 830nm can 
be efficiently detected with Silicon-based single photon 
counting modules. Even though KDP has a smaller non- 
linearity than other common x' 2 ' materials such as BBO, 
the long crystal used translates into a gain in production 
rate of photon pairs with respect to typical PDC sources 



As was discussed above, broadening of the joint spec- 
tral amplitude results from GVD terms, which [see 
Eqns. [25] are proportional to the crystal length. Thus, 
increasing the crystal length used leads to a temporally- 
broader two-photon time of emission distribution. If the 
mixed term phase 2/3 t + 4p in Eq. [55] does not van- 
ish, such broadening eliminates the "alignment" with the 
time axes. However, in this section we will show that, 
for high-aspect ratio spectrally factorizable two-photon 
states, in the long-crystal regime the mixed term respon- 
sible for introducing correlations in the times of emission, 
can be made arbitrarily small (even temporal broadening 
still occurs). In fact, an increase in crystal length gives 
a net reduction of temporal correlations. 

Let us assume that the crystal length and pump band- 
width are chosen so that the first condition for factoriz- 
ability (Eq. |2"9"]) is satisfied. Under these circumstances 
(see Eq. I31[) . the joint spectral amplitude may be ex- 
pressed as: 



exp 







\ V} 1 




exp 


(ra s ) 2 _ 







xexp [i{(3 s v 2 s + faff + P p v s v t )} (39) 



where we have expressed the idler width in terms of the 
aspect ratio r as <jj = ra s and where the GVD terms 
{fis, Pi and (3 P ) were denned in Eq. [25] Note that be- 
cause in the above expression the lack of factorizability 
resides in the phase term, the joint spectral intensity is 
factorizable. Nevertheless the presence of a mixed phase 
term, proportional to /3 P , translates into lack of factor- 
izability in the temporal domain. In order to study the 
temporal properties of the wavepacket, let us carry out a 
Fourier transform of the above joint spectral amplitude, 
to obtain the following joint temporal intensity: 



S(t s ,ti) = \T{t s ,U 



exp 



2tf 
5t 2 



exp 



2*f 

'at? 



exp [-2<T 2 M t s t(}0) 



where the temporal wavepacket width along t s is given 
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in terms of the reciprocal aspect ratio s 



"'by: 



St s = 



2v^/2 + <74(2/32 + /32 s 2 + 2/ a 



V) 



the temporal wavepacket width along U is given by: 



(41) 



su 



2V2J2 + o*(2ft + 0y + 20? s*) 



.S77. 



(42) 



and 



<T"S 



' M 



2(2 + ^(2/3| + %s*+2Pi a *)) 



(43) 



represents the mixed term coefficient. Specializing to a 
high aspect ratio spectrally factorizable state such as one 
produced by the long KDP cyrstal source discussed in 
the last section, the temporal walkoff between pump and 
the ordinary-ray photon r vanishes. Thus, according to 
Eq. [321 the ordinary-ray spectral width is then given by 
<7 S = tr where a represents the pump bandwidth. Like- 
wise, in the long-crystal limit, according to Eq. [32] the 
idler spectral width is given by l/(y/jT e ). This leads to 
the following expression for the reciprocal aspect ratio: 



(44) 



k' p )L 



Note, from the above expression, that the reciprocal as- 
pect ratio s is inversely proportional to the crystal length 
L. This fact, together with the linear dependence of each 
of the GVD terms (f3 s ,Pi and /3 P ) on the crystal length L 
can be used to express the mixed term coefficient cr^ in 
terms of its overall length dependence. We thus obtain: 



2 



4(1 + ^) 



a A k' p (k' s k' p ) 



^{k[-k' p f [i + a^(k' s -k' p y 



(45) 



The net crystal length dependence of a 2 M in the long- 
crystal regime is inverse quadratic. Hence increasing the 
crystal length does not lead to an increase in magnitude 
of the mixed term, but actually reduces the correlations. 
Recall that for a factorizable high aspect ratio two photon 
state, the pump bandwidth a and crystal length L do not 
have to be specified jointly i.e. there is no relationship 
that should be fulfilled between them to ensure factor- 
izability. Together with the dispersion insensitivity, the 
latter means that the crystal can be made arbitrarily long 
without introducing temporal correlations. The main re- 
maining considerations using a long crystal are excessive 
temporal broadening and absorption of PDC photons in 
the crystal. We note that if the pump field is chirped, 
(3 P should be substituted by 2(3 t + 2 m the expressions 



for the mixed term, Eqns. l43l and l45l (where fit is intro- 
duced in Eq. [25] represents the quadratic phase carried 
by the pump field) . The presence of chirp in the pump 
pulse can be treated similarly: upon increasing the crys- 
tal length L, the magnitude of the mixed term decreases. 
Thus, we conclude that (within the regime of validity 
of the approximations used) the source presented here is 
insensitive to those dispersive effects leading to the ap- 
pearance of temporal correlations, either introduced in 
the PDC crystal itself or indeed prior to the crystal. 

We can understand the dispersion insensitivity effect 
described above by realizing that dispersion effects in 
general require not only the presence of a dispersive phase 
but also broadband light. Engineering the joint spectral 
intensity by making it nearly monochromatic along ei- 
ther the signal or the idler axis, means that the result- 
ing narrowband photon will experience limited disper- 
sion. Furthermore, the mixed term (leading to temporal 
correlations) is also reduced by the monochromaticity of 
one of the photons. Fig. IIII Cl shows for the two-photon 
state produced by collinear degenerate type-II PDC in 
a 2cm-long KDP crystal cut at a phasematching angle 
of 68° yielding photon pairs centered at 830nm. The 
joint spectral intensity is shown in Fig. IIII Cf A) while 
the joint temporal intensity in Fig. IIII CT B). Because the 
parameters of the source fulfil Eq. [Ml the joint spectral 
intensity is factorizable. Furthermore, the joint tempo- 
ral intensity is also factorizable, for arbitrarily long crys- 
tals. Thus, the two-photon state produced by such a 
two-photon source is free from spectral (temporal) cor- 
relations, making it an ideal source of conditionally pre- 
pared pure and Fourier-transform limited single photon 
wavepackets. 
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FIG. 2: High aspect ratio two-photon state obtained by de- 
generate collinear PDC from a 2cm-long KDP crystal exhibit- 
ing space time decorrelation (A) Factorable joint spectral in- 
tensity (B) Factorable joint temporal intensity. 



IV. TWO PHOTON STATE ENGINEERING VIA 
CRYSTAL SEQUENCES WITH BI-REFRINGENT 
COMPENSATORS 

As has been shown in this paper, group velocity match- 
ing is a powerful technique which enables the engineer- 
ing of two photon states with desirable properties, tai- 
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lored for specific applications. Nevertheless, group veloc- 
ity matching crucially depends on the dispersion exhib- 
ited by the pump, signal and idler fields in the nonlinear 
crystal used, and tends to occur only in specific wave- 
length ranges. In this section we propose and analyse 
a technique in which the group velocity mismatch can 
be controlled at arbitrary wavelenghts. This comes at 
the cost, however, of an increased source complexity. In 
this scheme, a sequence of nonlinear crystals is used in- 
terspersed with birefringent spacers exhibiting a disper- 
sion such that the group velocity mismatch introduced 
by the crystal is compensated for by that in the spacer. 
This scheme can be thought of as analogous to band gap 
engineering using photonic crystals, except that group 
velocity, rather than phase velocity, is compensated. 




FIG. 3: This figure shows a schematic of the proposed crys- 
tal sequence with intermediate bi-refringent spacers. Each 
crystal has a length L while each spacer has a length h. 



phasematching function of a single crystal and the sec- 
ond factor incorporates the combined effect of the crystal 
and spacer dispersion. In order to carry out more explicit 
calculations, it is helpful to write down the crystal and 
spacer phasemismatch as a Taylor expansion [similar to 
that in Eq. 1231 however here we omit all terms of orders 
higher than 0(v 2 ) ]. Thus, we obtain, in terms of the 
frequency detunings = uj^ — luq (with fj, = s,i): 



$ = LAfc (0) + hAK {0) + T s u s + 



(50) 



where LAk^ and IiAk,^ denote the constant terms of 
the Taylor expansions for the crystal and spacer phase 
mismatch terms and 



[k' p (2cj ) - fc> )] L+ [«;(2wo) - <M] h (51) 



with fj, = s,i constitutes the first-order coefficients of 
the expansion. The term Afc(°) vanishes under the as- 
sumption that each of the crystals is aligned such that 
phasematching is attained. 

We can now analyze the influence of the proposed crys- 
tal and birefringent compensator assembly on the joint 
spectral correlations and introduce a new function 



Consider the experimental arrangement shown in 
Fig. [3] consisting of N identical crystals and N — 1 
linear £W spacers. Each of the crystals is assumed to 
be cut and oriented for degenerate collinear type-II PDC 
while the compensators are assumed not to exhibit a \^ 
nonlinearity. It is further assumed that each crystal has 
length L while each spacer has length h. The phase mis- 
match in each of the crystals is given by 



where k,. 



(46) 

denotes the wavenumber for 



^ (with n = p,s, i) 
each of the three fields taking into account dispersion in 
the crystal. The phase mismatch introduced by each of 
the spacers is equivalently given by: 

Ak — K p — K s — Ki. (47) 

where (with fi = p, s, i) now represents the wavenum- 
ber for each of the three fields taking into account dis- 
persion in the birefringent spacer. For an assembly of 
N crystals and N — 1 spacers the overall phasematching 
function can then be calculated as 



0/v 



(Ak, Ak) 



N-l 

E 

m=0 



e im(LAfc+hAK) sinc 



zil±sin(^) 



2_ 

sin(f ) 



-Ak 
2 

LAK 



sine 



where we defined the quantity $ as 
$ = LAk + hAn. 



(48) 



(49) 



Hence, apart from an overall phase factor the crystal as- 
sembly phasematching function is composed of the prod- 
uct of two distinct functions: one corresponds to the 



T N (x) 



1 sin(iVa;) 
N sin(a;) ' 



(52) 



which describes the modifications to the spectral struc- 
ture of the phasematching function [see Eq. 05] . As de- 
picted in Fig. IIVI the function Yn(x) exhibits for large 
values of N a periodic structure of narrow peaks, which 
are separated by 27r. While for odd values of N, all peaks 
are positive, for even values of N the peaks alternate 
between positive and negative values. The width of the 
peaks diminishes with increasing N. We may, thus, write 
the resulting phasematching function (where the phase 
term in Eq. [48] is to be neglected) for the crystal assem- 
bly as 



(53) 



where 4>(v s , v^) represents the single- crystal phasematch- 
ing function and xi^s, vf, N) the crystal assembly contri- 
bution. If we express the phasemismatch as Taylor series 
considering only terms up to the first order, we obtain 
for the phase matching function the contributions 



/iAk (0) + T s v s + T lV% 



and 



(i/ g , Vi) = sine 



TiVi) 



(54) 



(55) 



(see 

Eq. l5"Tj) , except for the absence of the spacer contribution: 



where the r M values are similar to [with fj, = s, i] 



[k' p {2Lu c )-k'^ c )]L. 



(56) 
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FIG. 4: Plots of the Tn(x) function showing different behav- 
ior for even and odd N; while for even N, peaks alternate with 
valleys (separated by 7r) for odd N, there are no valleys. The 
peaks and valleys become narrower with increasing N. Such a 
function describes the phasematching contribution from the 
crystal-birefringent spacer sequence. 



Let us now study the increased flexibility introduced 
into the two-photon state by such an assembly contribu- 
tion to the phasematching function. We start our analy- 
sis by noting that the contours of the x(^s, t^; N) function 
are straight lines oriented with a slope given by: 

tan = ~T i = ~(k> p -k> s )L + ( K > p - K ' a )h (57) 

Unlike for the case of a single crystal, where a change 
in crystal thickness merely results in an alteration of the 
spectral width (and further modifications to the phase- 
matching function cannot be attained except using differ- 
ent spectral regions or using different crystal materials), 
in the case of our crystal sequence source, the orientation 
of the spectral phasematching structure can be adjusted 
at will by varying the crystal to spacer thicknesses ra- 
tio L/h. The latter represents a very important added 
flexibility as it means that in principle an arbitrary ori- 
entation for the phasematching function orientation can 
be obtained. 

To illustrate the power of such engineering of PDC we 
will discuss in the following the specific case of synthesiz- 
ing a phasematching function with unit slope contours. 
As seen in section UlI Al (see also [TTJ| ) such unit slope can 
result in spectrally decorrelated two-photon states. For 
the crystal assembly the inverse group velocities k' , n'^ 



(with fi = evaluated at u>o and k' p , k' p evaluated at 
2wq, and the lengths L and h should satisfy the condition 



(fc; + fcj-2fcpi + («; + /<-2«;)ft 







(58) 



Eq.[58]closely resembles the group velocity matching con- 
dition [see Eq. [37] ■ I n fact, this condition tells us that 
unit-slope phasematching contours require that a gener- 
alized group velocity condition be fulfilled: the weighted 
average (with crystal and spacer lengths as weighting fac- 
tors) of the group velocity mismatch in the crystal and 
the spacer must vanish. Note that for this condition to 
be attainable, the spacer must exhibit a group velocity 
mismatch with the opposite sign as that of the crystal 
itself. A second condition which must be fulfilled relates 
to the constant phase term hAn^ [see Eq. [50] . If such 
a term differs from a multiple of 27r, the maxima of the 
Tjv(ic) function can shift away from degeneracy (indeed 
by tuning such constant phase term, the resulting modes 
shift on the v s — z/$ plane). Thus, a second condition on 
the spacer length which must be met is as follows: 



27rm 
Ak(°) 



(59) 



with m being an integer number. This tells us that the 
spacer length should be given as an integer multiple of the 
quantity 2tt/A«;( ), which limits the minimum allowed 
spacer thickness. Thus, if the conditions expressed by 
Eq ns . [58l and [59l are fulfilled, the crystal assembly phase- 
matching function can be written as: 



Xn(v s ,Vi) = T N (T-[v s - i>i]) 



(60) 



where: 



T_ = - ([#(wo) ~ k>(co )}L + [<M - 4(uo)]h) . 

(61) 

Such a function consists of "ridges" (alternated with 
"trenches" for even N) described by unit-slope contours 
with a FWHM (along the v s — Vi direction) given, in the 
limit of large N, by 



SX = 



ircNT- 



(62) 



with 72 = 1.39156.. (this numerical value results from 
calculating the width of the T n (x) function) and Ao de- 
notes the central PDC wavelength. The maxima of differ- 
ent "ridges" result from the argument of the Tjy yielding 
multiples of 2ir. In the case studied here where group 
velocity matching is attained [see Eq. [58], the spectral 
separation between subsequent ridges along the v — vi 
direction can be shown to be given by 



AA 



\2 
A 



%/2cT_ 



(63) 



Note that the generation of spectrally decorrelated state 
actually necessitates a source with a single "ridge" ; mul- 
tiple ridges yields a joint spectral intensity with multiple 
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peaks in v a — v± space which precludes spectral factoriz- 
ability. Because a single ridge is not attainable with 
our crystal assembly source, we aim to let the separa- 
tion between subsequent ridges be as as large as possible 
so that the portion of the quantum state exhibiting the 
desired behavior can be isolated, e.g. by weak spectral 
filtering. We would like to emphasize that such spectral 
separation between subsequent ridges Az/ is inversely pro- 
portional to the parameter T_, which in turn is propor- 
tional to the two thicknesses L and h. This leads to the 
manufacturing constraint that the individual crystal and 
spacer lengths should be as short as possible. As will be 
shown below this technique can indeed lead to a source of 
a decorrelated state with experimentally feasible crystal 
and spacer thicknesses. 

Altogether the resulting joint spectral amplitude for 
the two-photon state produced by such a crystal assembly 
is therfore given by: 



(64) 



where a{v s represents the pump envelope function, 
which is to be modelled as a Gaussian as in Eq. [201 Our 
approach is to let the bandwidth of the single crystal 
phasematching function (f>(v s ,Vi) be much larger than 
that of the function x(w s , oJi] N), resulting from the crys- 
tal assembly. Thus, the spectral structure imposed on the 
phasematching function by the crystal assembly domi- 
nates over that of a single crystal. Given that the band- 
width of individual ridges of the function xi^s, Vi\ n) de- 
pends inversely on N (while that of the function <p(v s , vi) 
does not depend on N), it is possible simply by using a 
large enough number of crystals to ensure that this con- 
dition on the widths of these two functions is fulfilled. 
Let us note that for the case of interest when the gen- 
eralized group velocity matching condition [see Eq. [58] 
is satisfied, the orientations of x( v s, Vf> N) and that of 
the pump envelope function are orthogonal to each other 
(along the v s + and v s — V{ directions respectively). 
This leads to the important conclusion that by choosing 
the pump bandwidth appropriately, a factorizable two- 
photon state can be obtained exhibiting nearly-circular 
contours. The value for the pump bandwidth a which 
guarantees such factorizable behavior can be expressed, 
in terms of the number of crystal segments N and the 
parameter T_ (see Eq. ET]) as: 
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In summary, the proposed strategy for synthesizing a 
decorrelated two-photon state is as follows: i) The crys- 
tal and spacer materials should be chosen such that the 
group velocity mismatch for the crystal and spacer mate- 
rial have opposite signs, ii) The crystal and spacer thick- 
nesses L and h should be chosen so that the generalized 
group velocity matching condition [see Eq. [55J is fulfilled 
and such that the constant phase term Ak' ' vanishes 
[see Eq. [SH] . Additionally, the crystal and spacer should 



be thin enough so that the spectral separation between 
modes of the xiys, v%\ N) function is large enough to en- 
sure that a single resulting mode may be isolated, iii) 
The number of crystal segments N should be made large 
enough so that the phasematching bandwidth of a single 
crystal is much larger than that of the assembly phase- 
matching function x(f s , Vi\ N). iv) The pump should be 
prepared so that it exhibits a bandwidth a which fulfils 
Eq. [65] thus guaranteeing that the pump envelope and 
phasematching functions are such that they yield a fac- 
torizable two-photon state. 




1.3 3,3 c . a 
signal wavelength l|±m 



signal wavelength lunnl 



FIG. 5: (A) Assembly contribution to the phasematching 
function for the case of 10 BBO crystals with a cut angle 
of 9 C = 42.35° (each with a thickness of 48.85mu) alternated 
with 9 calcite spacers (each with a thickness of 58.83/xm). 
Black areas indicate positive values while white areas indi- 
cate negative values. (B) Phasematching function for single 
BBO crystal, with above characteristics. (C) Pump envelope 
function for pump field centered at 400nm with a bandwidth 
of 1.48nm. (D) Resulting joint spectral amplitude function. 



We now study a specific example involving an assembly 
of BBO crystals and calcite spacers. Each of the crystals 
is characterized by a cut angle of 9 C = 42.35° yielding 
type-II collinear degenerate phasematching at a central 
PDC wavelength of A c = 800nm. We remark that BBO 
and calcite exhibit the required opposite sign group veloc- 
ity mismatch; for BBO: 2A£-A£-fc< = 3.535 x 10~ 4 ps/^m 
whereas for calcite: 2k' p — k'g — k^ = —2.936 x 10 _4 ps/ /im. 
The generalized group velocity matching condition [see 
Eq. [55] tells us that the ratio of the calcite spacer thick- 
ness to the BBO crystal thickness should be h/L = 1.204. 
This yields a value for the minimum spacer thickness 
27t/Ak(°) = 5.88/zm; we choose a value for the inte- 
ger to [see Eq. [52] of 10, thus yielding a spacer thick- 
ness h — 58.83/im and a BBO crystal thickness L = 
48.85/im. These thicknesses yield a spacing between as- 
sembly phasematching modes (along the v s — i>i direc- 
tion) of AA = 67.05nm, corresponding to a separation 
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along the v s or v± axes of 47.41nm. The latter means 
that over a bandwidth of 95nm around the central PDC 
wavelength (of 800nm) there is a single crystal-spacer as- 
sembly phasematching mode present. Fig. [5]JA) shows a 
plot of the function x(V s , Vi\ 10) for the above parame- 
ters, i.e. assuming that the assembly contains 10 BBO 
crystals and 9 calcite spacers. Note that near degeneracy 
(at A = 800nm), the slope of the contours defining the 
function is indeed unity. Fig. [DJB) shows a plot of the 
phasematching function 4>(u 8 , !/,•) for a single BBO crystal 
(with the above characteristics) making it clear that the 
spectral bandwidth exhibited is much larger than that 
of the assembly phasematching modes; the latter im- 
plies that the overall phasematching spectral structure is 
dominated by the assembly rather than the single crystal 
contribution. Fig. EJC) shows the pump envelope func- 
tion corresponding to a pump field with a center wave- 
length of 400nm and a spectral FWHM of 1.48nm, as re- 
quired by Eq. [55] to yield a factorizable state. Fig. [5]JD) 
shows the resulting joint spectral amplitude for our crys- 
tal assembly source. Due to the periodic character of the 
Tjv function there are multiple modes, separated along 
each frequency axis by 47 ' Anm. Nevertheless, weak spec- 
tral filtering can isolate the central mode (centered at 
As = Xi = 800nm). Fig. [6] shows a the joint spectral 
intensity corresponding to this central mode, clearly ex- 
hibiting spectral decorrelation. 

0.79 F 1 
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FIG. 6: This figure shows a the joint spectral intensity of the 
two-photon state produced by a sequence of 10 BBO crystals 
and 9 calcite spacers, as discussed in the text. The plot shows 
the region of the main degenerate peak near X a — Xi = 800nm. 

We have presented a novel technique which can yield 
factorizable states, based on group velocity matching 
compensation using a sequence of crystal segments alter- 
nated with bi-refringent spacers. The essential advantage 
of this approach is that it becomes in principle possible 
to synthesize a wide class of two-photon states simply 
by varying the spacer to crystal spacer thickness ratio, 
while choosing the crystal and spacer materials so that 
they obey certain relative dispersive properties. This 
scheme does not altogether eliminate the need for filter- 
ing; however, by reducing the crystal and spacer thick- 
nesses L and h (and thus increasing AA; see Eq. [55)1 
it becomes possible to increase the separation between 
modes, ideally to the point where in practice a single 



mode is present. In the example shown above, the thick- 
nesses were chosen in the region of 50/mi, since thinner 
segments may be difficult to manufacture and to handle. 
However, it is possible that by doping, ion exchange or 
some other process a single crystal may be modified in 
certain regions in such a manner that the group velocity 
mismatch changes sign; the modified regions would then 
play the role of the spacers, however in a single mono- 
lithic structure. With this approach it may be possible to 
take advantage of the smallest thicknesses allowed (with 
a value of m = 1 in Eq. |59"| . and thus eliminate the chal- 
lenges likely to be faced in assembling a large number of 
thin crystals. 



V. CONCLUSIONS 

We have studied single photon conditional preparation 
based on photon pairs generated by the process of para- 
metric downconversion. By modelling the detection as 
an appropriately defined projection operator we showed 
that, in general, the quantum state of the prepared state 
is described as a statistical mixture of the modes de- 
fined by the Schmidt decomposition of the PDC two- 
photon state. Furthermore, it was determined that a 
two photon state lacking correlations between the sig- 
nal and idler photons, in all degrees of freedom, is a 
basic requirement for ideal single photon preparation, 
i.e. yielding quantum-mechanically pure single photons. 
Such pure single-photon wavepackets are of fundamental 
importance in experiments relying on quantum interfer- 
ence of single photons from distinct sources, as is the case 
in the recently proposed scheme for quantum computa- 
tion with linear optics (4j. 

Assuming that spatial correlations can be indepen- 
dently eliminated (for example by the use of waveguided 
PDC@), we have discussed two techniques resulting in 
spectrally factorizable two-photon states. Firstly, it was 
shown that if an additional condition on the two pho- 
ton state is fulfilled, the technique reported by Grice et 
al. can yield a state in which both the joint tem- 
poral intensity, and the joint spectral intensity are fac- 
torable. The latter is crucial to guarantee full factor- 
izability in two photon states designed for conditional 
preparation of single photons. Secondly, we have intro- 
duced a novel technique in which the group velocity mis- 
match between the pump and downconverted photons, 
responsible for the mixedness of the prepared single pho- 
tons, can be controlled by using a sequence of crystals 
alternated with birefringent spacers exhibiting a disper- 
sion which compensates that of the crystal in a specific 
manner. This represents a powerful technique in which 
spectrally decorrelated states, as well as a more general 
class of spectrally engineered two photon states, can be 
obtained simply by varying the relative thicknesses of 
the crystal and spacer used. These techniques described 
may provide useful tools for practical implementations 
of novel quantum-enhanced technologies, such as linear 
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